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ABSTRACT 



Context. In this paper we analyze data of 8 elliptical galaxies in order to study the difference between their globular cluster systems 
(GCSs) radial distributions and those of the galactic stellar component. In all the galaxies studied here the globular cluster system 
density profile is significantly flatter toward the galactic centre than that of stars. 

Aims. A flatter profile of the radial distribution of globular cluster system respect to that of the galactic stellar component is a difference 
which has an astrophysical relevance. A quantitative comparative analysis of the profiles may give light on both galaxy and globular 
cluster formation and evolution. If the difference is due to erosion of the globular cluster system, the missing GCs in the galactic 
central region may have actually merged around the galactic centre and formed, or at least increased in mass, the galactic nucleus. An 
observational support to this is the clear correlation between the galaxy integrated magnitude and the number of globular clusters lost 
and that between the central massive black hole mass and the total mass of globular clusters lost. 

Methods. We fitted the stellar and globular cluster system radial profiles in a set of galaxies observed at high resolution. We saw that 
the globular cluster system profile is less peaked to the galactic centre than the stellar one is. Assuming this difference as due to GCS 
evolution starting from a radial distribution which was initially indistinguishable from that of stars, we may evaluate the number (and 
mass) of GCs 'disappeared' by a simple normalization procedure. 

Results. The number of missing globular clusters is significant, ranging from 21% to 71% of their initial population abundance in the 
eight galaxies examined. The corresponding mass lost to the central galactic region is 7 x 10 7 -1.85 x 10 9 M B . All this mass carried 
toward central galactic regions have likely had an important feedback on the innermost galactic region, including its violent transient 
activity (AGN) and local massive black hole formation and growth. 

Key words, galaxies: elliptical - clusters: globular 



1. Introduction 

Many elliptical galaxies contain more or less popolous globu- 
lar cluster systems (hereafter GCSs), that are, usually, less con- 
centrated towards the galactic centre than the bulge-halo stars. 
A huge amount of literature dedicated to GCS identification in 
external galaxies and study of their properties, since the sem- 
inal review by Harris & Racine (1979). Regard to the GCS ra- 
dial profiles, starting from the mentioned review, we just remind, 
for the sake of example, recent papers of Bassino et al. (2006), 
Goudfrooij et al. (2007), Lee et al. (2008), Rhode et al. (2007), 
Sikkema et al. (2006), Spitler, L.R. et al. (2006), Spitler, L.R. et 
al. (2008), Rhode & Zepf (2003), Peng et al. (2004) up to Harris 
et al. (2009). 

Though it cannot be yet concluded that this characteristic 
is common to all galaxies, we may say no case, at present, is 
known where the halo stars are less concentrated than the GCS. 
Moreover, there is a general agreement that the difference be- 
tween the two radial distributions is real and not caused by a 
selective bias. 

Consequently, different hypotheses have been advanced with 
the purpose of explaining this feature. Among these, two seem 
the most probable: 

(i) the difference between the two distributions reflects dif- 
ferent formation ages of the two systems, as suggested by Harris 
& Racine (1979) and Racine (1991); in their opinion globular 
clusters originated earlier, when the density distribution was less 



peaked. However, this hypothesis cannot explain why the two 
distributions are very similar in the outer galactic regions; 

(ii) another explanation is based on the simple assumption 
of coeval birth of globular clusters and halo stars, with a further 
evolution of the GCS radial distribution, while the collisionless 
halo stands almost unchanged. 

The GCS evolution is caused by both dynamical friction, 
which brings massive clusters very close to galactic centre, 
and tidal interaction with a compact nucleus (see for example 
Capuzzo-Dolcetta 1993). The combined effect of these dynami- 
cal mechanisms acts to deplete the GCSs in the central, denser, 
galactic regions, leaving the outer profile almost unaltered, then 
remaining similar to the profile of the halo stars. The efficiency 
of the mentioned phenomena is higher in galactic triaxial po- 
tentials (see for example Ostriker, Binney & Saha 1989; Pesce, 
Capuzzo-Dolcetta & Vietri 1992), in which there is a family of 
orbits, the 'box' orbits, which do not conserve any component 
of the angular momentum and then well populate the central 
galactic regions. Pesce et al. (1992) showed that globular clus- 
ters moving in box orbits lose their orbital energy at a rate an 
order of magnitude larger than those moving in loop orbits of 
comparable size and energy. These results showed the previous 
evaluations of the dynamical friction efficiency, based on very 
simplified hypotheses on the globular cluster orbit distribution 
that led to an overestimate of the dynamical braking time-scales. 
On the contrary, it has been ascertained that massive globular 
clusters in triaxial potentials are strongly braked during their 
motion and thus reach the inner regions in a relatively short 
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time (Capuzzo-Dolcetta 1993). There, they could have started 
a process of merging, giving origin to a central massive nucleus 
(eventually a black hole) or fed a pre-existent massive object 
(Ostriker et al. 1989; Capuzzo-Dolcetta 1993). Under the hy- 
pothesis that the initial GCS and halo-bulge radial distributions 
were the same, an accurate analysis of the observations would 
allow an estimate of the number of 'missing clusters' and there- 
fore of the mass removed from the GCSs. 

Actually, McLaughlin (1995), Capuzzo-Dolcetta & Vignola 
(1997), Capuzzo-Dolcetta & Tesseri (1999) and Capuzzo- 
Dolcetta & Donnarumma (2001), scaling the radial surface 
profiles of the halo stars of a galaxy to that of its GCS, estimated 
the number of missing globular clusters as the integral of the 
difference between the two radial profiles. Capuzzo-Dolcetta 
& Vignola (1997) and Capuzzo-Dolcetta & Tesseri (1999) 
suggested that the compact nuclei in our galaxy, M31 and M87, 
as well as those in many other galaxies, could have reasonably 
sucked in a lot of decayed globular clusters in the first few Gyrs 
of life. 

In this paper we enlarge the discussion of the comparison of 
GCS and star distribution in galaxies, dealing with eight galaxies 
for which good photometric data are available in the literature 
such to draw reliable radial profiles. 

In Sect. [2] we resume the way, discussed deeply in previously 
mentioned papers, to get from the compared GCS-stars radial 
profiles in a galaxy the number and mass of GCs disappeared due 
to time evolution of the GCS; in Sect. [3] we present and discuss 
the observational data, as well as the analytical fit expressions 
to the density profiles; in Sect. [4] we present the extension to the 
data of this work the correlation found by Capuzzo-Doletta & 
Donnarumma (2001) between the mass lost by the GCS, the host 
galaxy luminosity and the mass of the galactic central massive 
black hole. Finally, in Sect. [5] we summarize results and draw 
general conclusions. An error analysis of the methods used is 
presented in Appendix. 

2. The estimate of the number and mass of globular 
clusters lost 

Under the hypothesis that the flattening of the GCS distribution 
in the central region compared to the distribution of the stars 
in the galactic bulge, is due to an evolution of the GCS, the 
number of GC lost to the centre of the galaxy is obtained by 
the simple difference of the (normalized) density profiles inte- 
grated over the whole radial range (see McLaughlin 1995). A 
general discussion of the problems in evaluating the number and 
mass of missing clusters in this way can be found in Capuzzo- 
Dolcetta & Vignola (1997), Capuzzo-Dolcetta & Tesseri (1999) 
and Capuzzo-Dolcetta & Donnarumma (2001) and so it is not 
worth repeated here. 

In this paper, we assume as reliable fitting function for the GCS 
projected radial distribution a "modified core model", i.e. a law 



where So, r c and y are free parameters. The choice of this func- 
tion to fit the GCS profiles is motivated by the good agreement 
found with almost all the data used for the purposes of this pa- 
per, as confirmed by, both, the local maximum deviation of the 
fitting formula from the observed data and the computed^ 2 . For 
the galaxy stellar profile, S s (r), we rely on the fitting formulas 
provided by the authors of the various papers where we got the 



data from, that may change case by case, checking how good are 
the approximations to observed data. 

The "initial" distribution of GCs, EgcoMj (assumed to be equal 
in shape to the present stellar profile), can be practically ob- 
tained by a vertical translation of the stellar profile, E s (r), to the 
present GCS distribution, 2,cc( r ) as given by Eq.Q] We calculate 
the number of missing (lost) clusters as the surface integral of 
the difference between EgcoW and SccM over the radial range 
[0, r max ] of difference of these two profiles: 

Ni = 2n f "'" (E GC ,o(r) - S GC (r)) rdr. (2) 
Jo 

The present number of GCs, N, is obtained integrating Xcc( r ) 
over the radial range, [r m j n ,R], covered by the observations (for 
the value of R we rely on the papers where we got the GCs 
distribution data from). The values of N obtained with this 
method are usually different from those given by the authors of 
the papers, but for the purposes of this paper what is important 
is the difference bewteen N and Ni and, so, it is crucial a 
homogeneous way to determine them. 

The initial number of GCs in a galaxy is, indeed, estimated as 
^ =Ni + N. 

The numerical values of Ni, N, and N are functions of the fitting 
parameters and of the integration limits. These dependences and 
their contribution to the errors on the final results are discussed 
in Appendix. 

An estimate of the mass removed from the GCS (Mi) can be 
obtained through the number of GCs lost, Ni, and the estimate 
of the mean mass of the missing globular clusters, (mi). A 
priori, the determination of (mi) needs the knowledge of the 
initial mass spectrum of the CGS, which has suffered of an 
evolutionary erosion. However, the most relevant evolutionary 
phenomena (tidal shocking and dynamical friction) act on 
opposite sides of the initial mass function, and so we expect that 
the mean value of the globular cluster mass has not changed 
very much in time (see Capuzzo-Dolcetta & Tesseri, 1997 
and Capuzzo-Dolcetta & Donnarumma 2001). Hence, we can 
assume the present mean value of the mass of globular clusters, 
(m), as a good reference value for (mi). 

For NGC 4374, NGC 4406 and NGC 4636 we calculated (m) 
using their GC luminosity functions (GCLFs) and assuming 
the same typical mass-to-light ratio of GCs in our Galaxy, i.e. 
(MIL) %e = 1-5 for NGC 4406 and NGC 4636 or (M/L) B>0 = 1.9 
|25[) in the case of NGC 4374. For NGC 4636 we used also the 
mass function that represents its present distribution of GCs 
(see Sect. O- For the remaining galaxies (NGC 1400, NGC 
1407, NGC 4472, NGC 3258 and NGC 3268) there is no better 
way to estimate the total mass of 'lost' GCs than adopting as a 
'fiducial' reference value for their mean mass, (mi), the value, 
(m M w) = 3.3 x 10 5 M Q , of the present mean GC mass in our 
Galaxy. 



3. Data and results 

The data for the study of this paper have been collected from the 
literature. We will analyze a set of 8 galaxies for which the GC 
content and the radial profile has been reliably determined and 
apt to our purposes. 

The galaxies are: NGC 1400, NGC 1407, NGC 4472 (M 49), 
NGC 3268, NGC 3258, NGC 4374, NGC 4406 and NGC 4636. 
These galaxies go to enlarge the set of seventeen galaxies (Milky 
Way, M 31, M 87, NGC 1379, NGC 1399, NGC 1404, NGC 
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Fig. 1. Surface number density for NGC 1400, NGC 1407, NGC 4472 (M 49) and NGC 3268. Black squares represent the observed 
GC distribution; the solid line is its modified core model fit. The dashed curve is the surface brightness profile of the underlying 
galaxy (a power law and a central flat core for NGC 1400 and NGC 1407, a Sersic core model for NGC 4472 and a Nuker law for 
NGC 3268), vertically normalized to match the radial profile of the GCS in the outer regions. 



1427, NGC 1439, NGC 1700, NGC 4365, NGC 4494, NGC 
4589, NGC 5322, NGC 5813, NGC 5982, NGC 7626, IC 1459) 
whose GCS radial profiles have been compared to the stellar dis- 
tribution in previous papers (( 30), Capuzzo-Dolcetta & Vignola 
1997, Capuzzo-Dolcetta & Tesseri 1999, Capuzzo-Dolcetta & 
Donnarumma 2001). 



3.1. NGC 1400 

The surface density profile of this galaxy is given by Forbes et 
al. (2006) (hereafter F06) who fitted it by mean of a power law 
Y, s (r) oc r - 188 . This fitting law is reliable outside the galactic 
core, i.e. for r > rt = 0.0055 arcmin (Spolaor et al. 2008). The 



luminosity profile of the galaxy for r < r\, is almost flat, and 
linked to the external power law. We fitted the GCS distribution 
by a modified core model with So = 14.1 arcmin -2 , r c = 0.7 ar- 
cmin and y — 0.88. Integrating "Lcc( r ) m the radial range where 
GCs are observed, i.e. from r mm = arcmin to R = 2.8 arcmin, 
we obtain = 73 as the present number of GC in NGC 1400. 
The initial GCS distribution results to be approximated by: 

„ , . [ 1.37 x 10 5 arcmin -2 r < r b „. 

W) = ( 7 76r -1.88 r>rb 0) 

(see Fig.[T). 

Using the general method described in Sec. [2] and the estimated 
value r max - 2.3 arcmin we have that the number of missing 
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Fig. 2. Surface number density for NGC 3258, NGC 4374, NGC 4406 and NGC 4636. Black squares represent the observed GC 
distributions; the solid lines are their modified core model fit. The dashed curves are the surface brightness profile of the underlying 
galaxy (a Nuker law for NGC 3268, a power law and a central flat core for NGC 4374 and NGC 4636 and a Sersic core model for 
NGC 4406), vertically normalized to match the radial profile of the cluster system in the outer regions. 



clusters in this galaxy is Ni = 183, i.e. about 71% of the initial 
population of globular clusters, N, ■ = N/ + N = 256. 
An estimate of the mass lost by the GCS is M; = Ni (niMw) - 
6.04 x 10 7 M G . 



3.2. NGC 1407 

As for NGC 1400, data for this galaxy and its GCS are taken 
from F06. The luminosity profile of the galaxy stars is fitted by a 
power law, 2 s (r) oc r~ L42 . This power law fit fails in the inner re- 
gion where luminosity shows a core of radius rt, - 0.045 arcmin 
(Spolaor et al. 2008). As for NGC 1400 we thus assume, for 



r < r/,, a flat distribution matched to the external power law. The 
GCS modified core model has, in this case, So = 17.8 arcmin -2 , 
r c = 1.02 arcmin and y = 0.85 as better fitting parameters. The 
normalizing vertical translation of the stellar profile leads to 



v i \ _ \ 104 x 1° 3 arcmin 2 r < 

SgcoW - | 12 6r -1.42 



r> r b 



(4) 



(see Fig. [T}. 

Integrating the difference of the GCS "initial" and present radial 
profiles in the galactic region where these differ, i.e. up to r max - 
2.34 arcmin (see Eq. [2j, we obtain Ni = 84. The present number 
of GC obtained integrating I.cc( r ) in tne radial range covered by 
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the observations (i.e. from to R = 7.3 arcmin) is N = 314. The 
GCS has therefore lost 21% of its initial population, Ni = 398. 
Also in this case, we can evaluate the mass lost by the system as 
M, = N (m MW ) = 2.77 x 1O 7 M . 



3.3. NGC 4472 (M 49) 

Data for the GC distribution in this giant elliptical galaxy in 
Virgo are taken from Cote et al. (2003). The galaxy star lumi- 
nosity distribution, according to Ferrarese et al (2006), is well 
reproduced by a Sersic core model (Trujillo et al. 2004), i.e. by 

E,(r) = 74, {(7) V* - r) + e b "(^ 6(r - r h )eMiY J , (5) 

where 7 s (fb) = 7b, 0(x) is the usual Heaviside function, r b {break 
radius) divides the profile into an inner (r < r&) power-law region 
and an outer (r > r/,) exponential region; r e is the 'effective' 
radius and b n = 1.992« - 0.3271, with n free fitting parameter. 
For M 49 the parameter values are (Cote et al. 2003): y = 0.086, 
n = 5.503, b„ = 10.635, r e = 208 arcsec (= 16.89 kpc), r b = 
1.94arcsec (= 0.158 kpc). 

In this case the core model fit to the observed GCS distribution 
gives So = 6 kpc -2 , r c = 3.57 kpc and y = 0.652 as optimal 
parameters. The usual vertical translation leads to 

(/0.158\ 0086 
E GC , (r)=722j — — 0(0.158 -r)+ 



+94.1 G(r - 0.158)e" 10 - 635 ^)° 



(6) 



Fig.Q]shows how the modified core model profile does not repre- 
sent very well the observed GCS external profile, leading, there, 
to overestimate. 

The surface integral (Eq. [2]), performed with r max = 42.34 kpc, 
gives Ni = 5598. Integrating 7cc(r) up to R = 100 kpc we have 
that N = 6334. Hence, M 49 has lost 47% of the initial popula- 
tion of its GCs, Ni = 1 1, 932. 

As for NGC 1400 and NGC 1407, no better estimate of (mi) is 
available, so we evaluate M; as M; = N (m M w) = 1.85 x 1O 9 M . 

3.4. NGC 3268 

The GC distribution of this galaxy is discussed by Dirsch, 
Richtler & Bassino (2003). The resulting core model fit parame- 
ters have the values: Eo = 24.4 arcmin -2 , r c = 2.6 arcmin and 
y = 1.9. The stellar luminosity profile, following Capetti & 
Balmaverde (2006), is well represented by a "Nuker" law (in- 
troduced by Lauer et al. 1995): 



Z J (r) = 2^^(^) 7 



1 + 1^ 

rb 



(7) 



where B is the slope of the external region of luminosity profile, 
ri, is the "break radius" (corresponding to a brightness E/,) where 
the profile flattens to a smaller slope, measured by the parameter 
y; a sets the sharpness of the transition between the inner and 
outer profile. In the case of NGC 3268: a = 2.49, B = 1.64, y = 
0.13 and rb = 0.0252 arcmin. By the usual vertical translation of 
this profile, the initial GCS profile is obtained 



/0.0252V 

(— ) 



0.13 



1 + 



/ r \ 2 - 49 
1 0.0252/ 



-0.606 



The departure, visible in Fig.Q] of the profile given by Eqj8]from 
the modified core model profile in the external galactic region 
is mainly due to the incompleteness of the GC detection in the 
outermost regions. Integrating 7cc(r) from r m ,„ = arcmin to 
R = 7.94 arcmin we have N = 505 as present GC number. 
The number of globular clusters lost is found to be Ni = 398 
(from to r max = 2.3 arcmin), that means about 44% of the ini- 
tial abundance, Ni =913. 

Also in this case, to evaluate the mean mass of lost GC in 
NCG 3268 we have to assume (mi) = (niMw), obtaining M; = 
Ni(m MW ) = 1.31 x 1O 8 M . 

3.5. NGC 3258 

As for NGC 3268 the GCS density profile data for NGC 3258 are 
taken from Dirsch, Richtler & Bassino (2003). The best modi- 
fied core model fit is given by the values Eo = 16.9 arcmin -2 , 
r c = 3.1 arcmin and y = 2.4. The analytical fit to the luminos- 
ity profile of the galaxy is, again, obtained with the "Nuker" law 
(Eq. 0) with a = 2.10,/? = 1.51, y = and r b = 0.0192 arcmin. 
By mean of the usual procedure, the initial GCS profile is ob- 
tained: 

n-0.719 



EgcoM = 8489 



1 + 



/ r \ 210 
V0.0192/ 



(9) 



The present number of GCs is N - 343 (with R = 7.94 arcmin). 
Performing the surface integral of the difference of the initial and 
present distribution in the radial range up to r„ mx = 2.16 arcmin 
(see Eq.[2jl we have A^; = 212, corresponding to 38% of the initial 
GCS population, Ni = 555. 

For this galaxy, we obtained Mi = Ni (m M w) = 7.00 x 1O 7 M . 
3.6. NGC 4374 (M 84) 

Gomez & Richtler (2004) studied the GCS of this giant ellip- 
tical galaxy, using photometry in the B and R bands, to draw 
its radial surface distribution. Also in this case, the profile of 
the GC number density is flatter than the galaxy light (see Fig. 

The best modified core model fit to the GC data is given by 
Eo = 58.4 arcmin -2 , r c = 0.31 arcmin and y = 0.278. 
The galaxy light is characterized by a central core of radius 
rb — 0.0398 arcmin (Lauer et al., 2007); for r > rb, it is well 
fitted by the power law E. s (r) cc r -1 67 (Gomez & Richtler 2004). 
The usual normalization leads to 



Ecco(r) 



2.94 x 10 4 arcmin 2 r < r b 
135r -L67 r > r b , 



(10) 



(8) 



as GCS initial radial profile. Integrating our core model up to 
R - 11.8 arcmin we get Af = 4655 as present number of 
GCs. The usual integration of the difference of the initial and 
present GC distribution (Eq.|2]with r max = 3.84 arcmin) leads to 
N = 2361. Hence NGC 4374 has lost 34% of its initial popula- 
tion of globular clusters, Af, = 7016. 

In the case of NGC 4374 the value of the mean mass of a GC has 
been evaluated using the GCLF in the R band given by Gomez & 
Richtler (2004). The mean color ((B - R) ) = 1 . 1 8 of GCs in this 
galaxy (Gomez & Richtler 2004) allows us to estimate the mean 
B absolute magnitude and the mean luminosity of GCs in the B 
band, ((L/Lb) ) assuming m - M = 31.61 (Gomez & Richtler 
2004). It results ({L/L B ) Q ) = 1.75 x 10 5 (withM B0 = 5.47, Cox 
2000). 

Adopting the mass to light ratio (M/L)b.q) - 1.9 obtained 
by Illingworth (1976) for 10 galactic globular clusters, we get 
(mi) = 3.33 x 1O 5 M and M, = N t (m t ) = 7.86 x 10 8 M o . 
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Galaxy £ r c y r max R 



NGC 1400 


14.1 


0.7 


0.88 


2.3 


2.8 




1 7 S 
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i m 
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U.oJ 




7 ^ 

1.3 


NGC 4472 


142 


0.73 


0.652 


8.69 


20.52 


NGC 3268 


24.4 


2.6 


1.9 


2.3 


7.94 


NGC 3258 


16.9 


3.1 


2.4 


2.16 


7.94 


NGC 4374 


58.4 


0.31 


0.278 


3.84 


11.8 


NGC 4406 


26.76 


3.52 


1.19 


5.36 


24 


NGC 4636 


77.66 


0.823 


0.691 


4.75 


6.6 



Table 1. Col. 1: galaxy name; col. 2, 3 and 4: parameters of the modified core model fit for all the galaxies studied; col. 5: upper 
limit in the integral giving the number of the lost GCs (r max ); col. 6: upper limit of the integral performed to estimate the present 
number of GCs (R). So is in arcmin" 2 ; r c , r max and R are in arcmin. 



Galaxy 


Model 


1 






b„ 


n 


a 




J 


NGC 1400 


lc 


7.76 


5.5 x 10" J 








1.88 






NGC 1407 


lc 


12.6 


0.045 








1.42 






NGC 4472 


cS 


17139 


0.0323 


3.47 


10.635 


5.503 






0.086 


NGC 3268 


N 


13692 


0.0252 








2.49 


1.64 


0.13 


NGC 3258 


N 


8489 


0.0192 








2.10 


1.51 





NGC 4374 


lc 


135 


0.0398 








1.67 






NGC 4406 


cS 


13490 


0.012 


6.86 


13.649 


7.02 






0.021 


NGC 4636 


lc 


70.8 


0.0573 








1.5 







Table 2. Galactic luminosity fitting parameters. Col. 1 : galaxy name; col. (2) key identifying the galaxy light profile model (lc=linear 
with a flat core in the inner region, cS=core-Sersic, N=Nuker); col. 3-10: parameters of the various profile models (see Sect. [2] and 
[3]for details), rj is in arcmin" 2 ; and r e are in arcmin. 



3.7. NGC 4406 (VCC 881) 

NGC 4406 is another giant elliptical in Virgo; its GCS has been 
studied by mean of the Mosaic Imager on the 4m Mayall tele- 
scope at the KPNO (Rhode & Zepf 2004) in the B, V and 
R bands. The resulting best fit core model is characterized by 
So = 26.76 arcmin" 2 , r c = 3.52 arcmin and y — 1.19. The galaxy 
light is well fitted by a Sersic core model (Eq.[5]l, whose parame- 
ters have been determined by Ferrarese et al. (2006). Its vertical 
translation gives the initial GCS radial profile 



SccoM = 13490 



0.012 



0.021 



0(0.012 - f) 



+ 250(9(r - 0.012)e" 13 - 649 ^)° 



(11) 



Integrating the present distribution of GCs, from r mm = arcmin 
to R = 24 arcmin, we have N = 2850. The surface integral given 
in Eq.|2] with r max = 5.36 arcmin, gives the number of globular 
clusters lost, Ni = 1359, i.e. about 32% of the initial GC popula- 
tion. 

Using the GCLF of this galaxy and its distance modulus m-M = 
31.12 (Rhode & Zepf, 2004), we evaluated the mean value of the 
absolute GC V magnitude, (M v ) - -8.42 which corresponds to 
the mean luminosity (L/L Q ) V = 1.98 x 10 5 (My® = 4.82 from 
Cox 2000). 

Assuming (M/L)y >0 = 1.5, we obtain (m/) = 2.97 x 10 5 M o . 
This estimate leads to the value of the mass lost by the GCS, 
Mi = Ni (mi) = 4.04 x 1O 8 M . 

3.8. NGC 4636 

The GC content of this galaxy has been studied by Kissler Patig 
et al. (1994). The modified core model fit has So = 77.66 
arcmin" 2 , r c = 0.823 arcmin and y = 0.691 as optimal parame- 



ter values. 

The galactic light profile shows an inner flat distribution (a core 
with radius r\, ^ 0.0573 arcmin (Lauer et al. 2007)), while for 
r > r/, the light distribution is well fitted by the power law fit 
S,(r) oc r" 1 5 (Kissler Patig et al. 1994). 

The vertical translation of the stellar profile gives the initial GCS 
profile: 



£cc,o(r) = 



5.16 x 10 arcmin r < r b 
70.8r" 15 r > r b . 



(12) 



Integrating the present surface density profile of the GCS up to 
R = 6.6 arcmin, we obtain N = 141 1. Performing the surface in- 
tegral given in Eq.|2](with r max = 4.75 arcmin), we estimate that 
the number of GCs disappeared is N/ = 746, i.e. 35% of the ini- 
tial population, A^, = 2157. In the case of this galaxy we obtained 
two different estimates of the mass lost by the GCS, starting from 
data taken from Kissler Patig et al. (1994). The first estimate has 
been obtained using the GCLF (Kissler Patig et al. 1994). As 
for NGC 4406 we calculated the mean absolute V magnitude of 
GCs, (M v ) = -8.07, (given m - M = 3 1 .2 by Kissler Patig et al. 
1994). Assuming for GCs in NGC 4636 the same M/L v ratio of 
galactic GCs, (M/Ly)® = 1.5, the deduced mean luminosity of 
GCs, (L/L B )y = 1.43 x 10 s , gives (m u ) = 2.15 x 1O 5 M , and 
so M u = N, <m/,i) = 1.29 x 10 8 M o . 

Another estimate is found using the mass distribution of GCs 
obtained in Kissler Patig et al. (1994) transforming the magni- 
tude bins of the GCLF candidates into masses using the relation 
given by Mandushev et al. (1991): log(M/M ) = -0A6M V + 1.6 
(corresponding to a mean mass to light ratio (M/L)y ^ 2.0). 
Knowing the mass distribution we can directly calculate the 
mean mass of GCs, (m/ j2 ) = 3.79 x 1O 5 M , and thus M;, 2 = 
Ni(m L2 ) = 2.97 x 1O 8 M . 

The averages of our two estimates gives Mi = 2.22 x 1O 8 M . 
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Tables [T] presume the parameters of the radial profile fitting 
functions for the studied galaxies, while [3] resume the results in 
terms of estimated number and mass of GC lost. 

4. The correlation between M t , M v and M bh 

The evolutionary explanation of the difference between the ini- 
tial and present GC distribution implies a correlation between 
the (supposed) mass lost by GCS with the mass of the galactic 
central supermassive black hole (M/,/,) and. likely, with the host 
galaxy luminosity (My). Tab.|4]reports the whole set of galaxies 
for which we have the estimate of My, Mbh and M/. 
Fig. [3] shows a plot of Tab. [4] data which clearly indicate an in- 
creasing trend of M/ as function of Mbh (left panel) and of My 
(right panel). 

In particular, the linear fit of data in the left panel is given 
by log Mi - a log M b h + b with a ± <r(a) = 0.47 ± 0.20 and 
b + cr(b) = 3.9 + 1.8, giving r = 0.45 and x 2 = 9.5. The, alter- 
native, exponential fit gives log Mi - a exp(log M/,h) + f3 where 
a+o-(a) = (7.7±3.3)xlO- 5 andy6+0-OS) = 7.45+0.28, r 2 = 0.20 
and x 2 = 9.6. 

We performed also fits for for the correlation between M/ and 
Mbh excluding the data which have a great residual from the 
mean square (those of NGC 1439 and NGC 1700) obtaining the 
following parameters of the linear fit: a + cr(a) - 0.61+0.15 and 
b±cr(b) = 2.9 + 1.3, with r - 0.69 and* 2 = 4.0. The exponential 
fit parameters are in this case: a + cr(a) = (1.17 + 0.22) x 10~ 4 
and p + cr(J3) = 7.31+0.18, giving r 2 = 0.60 and X 2 = 3.1. 

The least square, straight-line fit to the whole set of data 
shown in the right panel of Fig. [3] is given by log M/ = aMy + b 
with a + cr(fl) = -0.62 + 0.15 and b + cr(b) = -5.3 + 3.2, giving 
r - 0.67 and* 2 = 6.6. The exponential fit on the same data gives 
log Mi = a exp(-M v )+j3 where a + cr(a) = (2.20 + 0.49)x 10" 10 
and [5 + cr(/3) = 7.33 + 0.19, r 2 = 0.48 and * 2 = 6.2. 

The correlation seen in the right panel of Fig. |3]between Mi 
and My reflects, both, an expected physical dependence on the 
total galactic mass of evolutionary processes acting on GCSs 
and, simply, the positive correlation between Mbh and My for 
the same set of galaxies. Actually, the MbhrMy correlation for 
the set of galaxies in Tab. |4] has a clearly positive slope, as 
shown also by the least square fit in Fig. [4] The least square 
fit is log M b h = aM v + b with a + <r(a) - -0.56 + 0.15 and 
b + cr(b) = -3.41 + 3.22, giving* 2 = 6.54. On the other 'phys- 
ical'side, the energy and angular momentum dissipation caused 
by dynamical friction should depend on the inner galaxy phase 
space density (oc p/cr 3 , where p and cr are the galactic mass den- 
sity and velocity dispersion, respectively). A stronger dynami- 
cal friction causes a faster GC decay toward inner galactic re- 
gions where the tidal action of a massive black hole depletes 
the GC population. Were brighter galaxies also denser in the 
phase-space, the Mbh ~ My correlation would have a p/cr 3 vs. 
My counterpart. Using data available in the literature for a set 
of 428 galaxies (the largest part coming from a combination of 
data available in the catalogue by Prugniel & Simien (1996)) we 
find the distribution shown in the right panel of Fig. [4] which, 
far from being conclusive, shows indeed a trend of higher cen- 
tral phase space density in brighter galaxies. The least square fit 
is log(p/cr 3 ) = aM v + b with a + aid) = -0.07 + 0.019 and 
b + oib) = -1.65 ± 0.39, giving* 2 = 194.44. 



eight elliptical galaxies observed by various authors. We find 
that GCS distributions flatten toward the centre, showing a broad 
core in the profile, contrarily to the surrounding star field. This 
result agrees with many previous findings, indicating, indeed, 
that GCs are usually less centrally concentrated than stars of the 
bulge-halo. A debate is still open on the interpretation of this 
observational issue. The "evolutionary" interpretation is partic- 
ularly appealing; it claims that, initially, the GCS and stellar 
profiles were similar and, later, GCS evolved to the presently 
flatter distribution due to dynamical friction and tidal interac- 
tions (Capuzzo-Dolcetta 1993). In this picture, the flatter cen- 
tral profile is due to the erosion of the inner GCS radial profile. 
Many GCs are, consequently, packed in the inner galactic re- 
gion, where they influence the physics of the host galaxy. Many 
of the galaxies studied so far have massive black holes at their 
centres, whose mass positively correlates with our estimates of 
number and mass of GC lost. This is a strong hint to the validity 
of the mentioned evolutionary scenario, together with the other 
evident correlation between number and mass of GC lost and 
their parent galaxy luminosity. The evolutionary hypothesis is 
also supported by the positive (although statistically uncertain) 
correlation between the (rough) estimate of the galactic central 
phase-space density and integrated magnitude. At the light of 
these encouraging findings, we think that much effort should be 
spent into deepening the observational tests of this astrophysical 
scenario. 



Appendix A: The error on the estimates of number 
of lost GCs 

Here we describe how we evaluated the errors, e/, given in Table 
[3] As explained in Sect. [2] the number of GCs lost in the galaxies 
of the sample has been evaluated as the integral of the difference 
between the (estimated) initial and present GCS radial distribu- 
tions over the radial range [r,„,„, r max ] where the two profiles dif- 
fer. The absolute errors on Ni (AN/) are given by the sum of the 
error on Ni (ANt) and the error on (AN), yielding the relative 
error e; = of Table [3] We may estimate AN and A^V, as fol- 
lows. 

i) estimate of AN 

For all the galaxies the number N(r min , r max ) is given by 



N(r m i n ,r max ) — 27r£i 



Z 7T(r 2 + r 2 ) 



(1-7) 



(A.l) 



which is a function of the parameters So, r c , y, r m ,„ and r„, 
whose indetermination is 



AN = 



dN 



<92o 
dN 





dN 




dN 


AE„ + 


Ar c + 


dr c 


dy 



Ay + 



dffnbi 

where: 





dN 


min 




max 



Ar„, 



(A.2) 



5. Conclusions 

We presented the comparative discussion of radial distribution 
of the globular cluster systems and of the stars in a sample of 



1 + 



(*)' 



-dr ■ 



n(r 2 + r?) 



(i-r) 



l + 



ar 



,(A.3) 



8 



R. Capuzzo-Dolcetta and A. Mastrobuono-Battisti: Globular Cluster System erosion in elliptical galaxies 





Fig. 3. The correlation between the GCS (logarithmic) mass lost and the central galactic black hole mass (left panel) and integrated 
V magnitude of the host galaxy (right panel) for the set of galaxies in Table |4] Masses are in solar masses. Black circles represent 
the eight galaxies whose GCS data are discussed in this paper, black triangles refer to the others. The straight lines and curves are 
the approximation fits discussed in Sect0] 



m _ r- 



27vY Q (r 2 c + yr 2 ) 



(l+y) 



dr ■■ 



dr — 



7rS () (r 2 + r c 2 ){l+( r -l)ln 




(r-D 2 




y 



dN 



— — 27rL 



(we set r m j n = 0.1 arcmin). 



dN 



dr 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



For NGC 1400, NGC 1407, NGC 4374, NGC 4636 we have (see 
Sect. 13. II 13.21 13.61 l3~8l and Tab.[2]for the meaning and the values 
of the parameters) 

Ni(r min ,r max ) = 2nr]r h a f rdr + 2nr] f r 1 ~ a dr = 



— 7ri]r b a r 2 \'^ + 2nrj 



Jl-a 



2-a 



(A.8) 



In Eq. IA.8I r\ represents the parameter obtained by the vertical 
shifting of the luminosity profile. 
The error AN, is thus given by: 



ANi 



dr] 
where: 



At] ■ 



8Ni 



da 



Aa ■ 



dN, 



dr„ 



At" min 



dN 



max 



Ar max (A.9) 



dN- C b ( Vmax 

— - = 2;rrr a rdr + 2n r x ~ a dr ■■ 
dV " X„„„ L 



-a ,2\<~b 



Jl-a 



+ 2n- 



2-a 



dN-, 
dr b 

dN 

The fitting parameters used to calculate AN are summarized in ^ a 
Tab. Q] 

ii) estimate of ANj 

The fitting formulas to the initial distribution of GCs change 
for the various galaxies studied. 



rdr = —nrjr^^r 2 ^ 



r> 

= —2nj]r b a I rln(ri,)dr - 2ttt] 
= -2nr, a r 2 \n(r h )\[\ + 



2m 1 r 2 - a {\ + [a -2] ln(r)} 



•Jr,, 



( r -2) 2 



(A.10) 
(A.ll) 

ln( r)dr = 

(A. 12) 
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Fig. 4. Left panel: the correlation between the logarithm of the central galactic black hole mass (in solar masses) and the integrated 
V galactic magnitude (see Table 0J. Right panel: the correlation between the value of the parameter proportional to the central 
galactic phase-space density (in arbitrary units) and the galactic integrated V magnitude. The straight lines are the least-square fits 
to the data (see Sect.g]) 



Also in this case we assumed r m ,„ =0.1 arcmin. For all the galax- dN, 
ies analyzed, r m ,„ > r b and r max > r b ; so we have 



dNi 

3? min 

dNj 

$F max 



= -2nr\r m *, 



db„ J n 
(A-13) dNi 



(A. 14) 



For M 49 and NGC 4406 (Sect.[33]and Sect. [3/7]) we have 

Nj(r m i n , ^max) ~ 



2nn 



6(r b -r) + e r e ' \0{r - r h ) ) dr ■■ 



2nrj 2 ( r b ' 



dr (A. 15) 



where b„ = 1.992n - 0.3271. The second row of the previous 
expression is justified by the fact that, both for M 49 and NGC 
4406, r b > r mm . Thus the error on N\ 



AN/ 



dNi 


An + 


dN 


Ab n + 


dN 


Ar h + 


dN 


dr/ 


db n 




dr b 




dr e 



Ar e + 



dNi 


An + 


dNi 


min 


dN/ 


dn 




3^ min 




$F max 



Ar„ 



(A. 16) 



where Ab„ = An, is evaluated by the following expressions 
of the individual error contribution: 



m 

dn 



2n 2 (r b Y 
-r 



2-7 



re 1 Jar 



+ 2,£\^W-m d , (A.17) 



re 1 



-f -\J 

r e \r e 



dr, (A. 18) 
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dr e nrj J n 
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dn 



-In 
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dr. 



(A.19) 
(A.20) 

(A.21) 



(A.22) 



(A.23) 
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Remembering that, for both NGC 4472 and NGC 4406, r„„„ > r b 
and r max > r b we can estimate the following contributions 



(A.24) 




(A.25) 



See Tab.|2]for the values of the parameters used in the Eq. IA.15l - 
Ea EA25l 

Last for NGC 3268 and NGC 3258 (see Sect. 03 Sect. [53] References 
and Tab. [2] for the values of the parameters) we have 



All the integrals from Eq. |A.15| to Eq. |A.32l must be calcu- 
lated numerically using the values of the parameters given in 
Table [2] The results listed in Table Q] are obtained assuming an 
error of 1 % on each independent parameter used. Only in the 
case of NGC 3258 we assumed Ay = 0.001 because Capetti & 
Balmaverde (2006) obtained y — from their Nuker fit, and so 
it is impossible to give an estimate of the error as a percentage 
of y. 
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0.38 


0.16 


1.83 x 10 8 


7.00 x 10 7 


NGC 4374 


4655 


7016 


2361 


0.34 


0.050 


2.34 x 10 9 


7.86 x 10 8 


NGC 4406 


2850 


4209 


1359 


0.32 


0.23 


1.25 x 10" 


4.04 x 10 8 


NGC 4636 


1411 


2157 


746 


0.35 


0.11 


6.41 x 10 8 


2.22 x 10 8 



Table 3. col. (1): galaxy name; col. 2-8: the present number of GCs (AO, its initial value (Ni), the number of GCs lost (Ni), the 
percentage of GCs lost and the estimated relative error on Ni (e;, see Appendix), the estimate of the initial mass of the whole GCS 
(Mi) and of the mass lost by each GCS (Mi). Mi and Mi are in solar masses. 



Galaxy 


M 


V 






M t 


Sources 


MW 


-20 


.60 


3 


.61 


X 


10" 


1.80 x 10' 


E05, WOO, CDV97 


M31 


-19 


.82 


6 


.19 


X 


10 7 


2.30 x 10 7 


M98 , CDV97 


M87 


-22 


.38 


3 


.61 


X 


10 9 


2.33 x 10 9 


M98, CDV97 


NGC 1427 


-20 


.43 


1 


.17 


X 


10 s 


8.86 x 10 7 


VM99, CDT99 


NGC 4365 


-22 


.06 


7 


.08 


X 


10 8 


7.48 x 10 7 


VM99, CDT99 


NGC 4494 


-20 


.94 


4 


.79 


X 


10 8 


2.98 x 10 7 


VM99, CDT99 


NGC 4589 


-21 


.14 


3 


.09 


X 


10 8 


7.58 x 10 7 


VM99, CDT99 


NGC 5322 


-21 


.90 


9 


.77 


X 


10 8 


-6.51 x 10 7 


VM99, CDT99 


NGC 5813 


-21 


.81 


2 


.82 


X 


10 8 


1.03 x 10 8 


VM99, CDT99 


NGC 5982 


-21 


.83 


7 


.94 


X 


10 8 


8.86 x 10 7 


VM99, CDT99 


NGC 7626 


-22 


.34 


1 


.95 


X 


10 9 


3.59 x 10 7 


VM99, CDT99 


IC 1459 


-21 


.68 


2 


.60 


X 


10" 


1.57 x 10 8 


Fe05, VM99, CDT99 


NGC 1439 


-20 


.51 


1 


.95 


X 


10 8 


4.79 x 10 6 


VM99, CDT99 


NGC 1700 


-21 


.65 


4 


.37 


X 


10 9 


3.66 x 10 6 


VM99, CDT99 


NGC 1399 


-21 


.71 


5 


.22 


X 


10" 


1.44 x 10 8 


M98, CDD01 


NGC 1400 


-20 


.63 


3 


.71 


X 


10 8 


8.45 x 10 7 


VM99, F06, CDM 


NGC 1407 


-21 


.86 


5 


.55 


X 


10 8 


1.31 x 10 8 


Z07, F06, CDM 


NGC 4472 


-23 


.10 


2 


.63 


X 


10" 


3.94 x 10 9 


M98, RZ04, CDM 


NGC 3268 


-22 


.07 


4 


.68 


X 


10 8 


2.98 x 10 8 


BC06, D03, CDM 


NGC 3258 


-21 


.40 


2 


.14 


X 


10 8 


1.83xl0 8 


BC06, D03, CDM 


NGC 4374 


-22 


,62 


1 


.41 


X 


10 9 


2.34 x 10 9 


R98, GR04, CDM 


NGC 4406 


-22 


.30 


1 


.40 


X 


10 8 


1.25 x 10 9 


CJ93, RZ04, CDM 


NGC 4636 


-21 


.70 


3 


.63 


X 


10 8 


6.41 x 10 8 


VM99, KR94, CDM 



Table 4. col. (1) galaxy name; col. (2), (3) and (4): the V absolute magnitudes, the galactic central black hole masses and the mass 
lost by GCSs (both in solar masses), respectively; col. (5): bibliographic reference sources for entries in col. (2), (3) and (4); CDM 
is the present paper, the other ackronyms are defined in the References. 



